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ABSTRACT 

We study the chemical and spectro-photometric evolution of galactic disks with 
detailed models calibrated on the Milky Way and using simple scaling relations, based 
on currently popular semi-analytic models of galaxy formation. We compare our re- 
sults to a large body of observational data on present day galactic disks, including: 
disk sizes and central surface brightness, TuUy-Fisher relations in various wavelength 
bands, colour-colour and colour-magnitude relations, gas fractions vs. magnitudes and 
colours, abundances vs. local and integrated properties, as well as spectra for different 
galactic rotational velocities. Despite the extremely simple nature of our models, we 
find satisfactory agreement with all those observables, provided the timescale for star 
formation in low mass disks is longer than for more massive ones. This assumption is 
apparently in contradiction with the standard picture of hierarchical cosmology. We 
find, however, that it is extremely successfuU in reproducing major features of present 
day disks, like the change in the slope of the TuUy-Fisher relation with wavelength, 
the fact that more massive galaxies are on average "redder" than low mass ones (a 
generic problem of standard hierarchical models) and the metallicity-luminosity rela- 
tion for spirals. It is concluded that, on a purely empirical basis, this new picture at 
least as successful as the standard one. Observations at high redshifts could help to 
distinguish between the two possibilities. 
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1 INTRODUCTION 

Considerable progress has been made in the past ten years 
or so towards a qualitative understanding of galaxy prop- 
erties. This has been made in the framework of the Cold 
Dark Matter (CDM) picture for galaxy formation, either 
with full N-body simulations coupled to gas hydrodynamics 
(e.g. Steinmetz 1998) or with simpler, so-called semi-analytic 
(SAM) models (e.g. Somerville and Primack 1998). 

N-body simulations have difficulties up to now in re- 
producing the observed properties of galaxies in detail. In 
particular they fail to produce realistic models for present 
day spirals, because excessive transfer of angular momentum 
from gas to dark matter results in overly small disks (e.g. 
Navarro and Steinmetz 1997). Moreover, the first attempt 
to combine in a single computation dynamical, chemical and 
photometric evolution (Contardo et al. 1998) failed to repro- 
duce basic observational properties of the Milky Way disk 
(hke the local G-dwarf metallicity distribution or colour gra- 
dients at low redshift). 

Semi-analytic models, pioneered by White and Frenk 
(1991) and subsequently augmented by Monte-Carlo meth- 
ods, were developed by several groups ("Munich": Kauff- 



mann. White and Guiderdoni 1993, Kauffmann 1996, Kauff- 
mann, Nusser and Steinmetz 1997; "Durham": Cole et al. 
1994, Baugh et al. 1997; "Santa-Cruz": Somerville and Pri- 
mack 1998). Such models include a simplified treatment of 
several essentially unknown ingredients, like star formation, 
gas cooling, supernova feedback and galaxy merging. They 
constitute an efficient tool to explore the large parameter 
space occupied by these (and several other!) unknowns, but 
they also provide an important level of understanding that 
could not be achieved by running N-body simulations alone. 
The successes and shortcomings of such models are nicely 
reviewed in Somerville and Primack (1998). 

Despite the apparent success of SAMs, it is not obvi- 
ous whether their number of free parameters is smaller than 
the number of the observables they successfully explain. One 
rarely mentionned shortcoming is that the assumption of In- 
stantaneous Recycling Approximation (IRA) results in over- 
estimated metallicities and underestimated gas fractions at 
late times. Moreover, the wealth of available spatial informa- 
tion (i.e. data on gradients of colours, gas and metallicity) 
as rarely used as constraints to SAMs (with the exception 
of Kaufmann 1996). 
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This latter point is dealt with in Jimenez et al. (1998), 
who use the simple scaling relations for galactic disks (es- 
tablished in the framework of SAMs, see Sec. 2.2) in order 
to derive initial gaseous profiles. They subsequently calcu- 
late the full chemical and spectrophotometric evolution of 
such disks for the case of Low Surface Brightness (LSB) 
galaxies and compare succesfuUy their resulting colour and 
abundance gradients to observations. 

In this paper we adopt an approach similar to that of 
Jimenez et al. (1998), i.e. we use multi-zone models and sim- 
ple scaling laws given by SAMs in order to study the evolu- 
tion of galactic disks. Our method and purpose are, however, 
different. We calibrate our models to a fairly detailed "tem- 
plate" of the Milky Way (elaborated in a previous paper: 
Boissier and Prantzos 1999, paper I in this series), while 
Jimenez et al. (1998) make no such calibration. Moreover, 
we use a different prescription for the radial dependence of 
the star formation rate, inspired by theories of star forma- 
tion in spiral disks (see Sec. 2.1). Finally, we are interested 
mainly on the properties of High Surface Brightness disks, 
although (as we shall see) some of our models produce LSBs 
and some bulges or ellipticals. 

Our main objective is to check whether in such an over- 
simplified framework (i.e. neglecting all the poorly known 
physics of star-gas interactions, galaxy merging etc.) one can 
reproduce observations concerning present day disks. For 
that purpose we use a large body of observational data on 
disks, concerning sizes and central surface brightness, TuUy- 
Fisher relations in various wavelength bands, colour-colour 
and colour-magnitude relations, gas fractions vs. magnitudes 
and colours, abundances vs. local and integrated properties, 
as well as spectra for different rotational velocities. Compari- 
son of the models to the data turns out to be successfull, pro- 
vided a crucial assumption is made: that small galaxies form 
on average their stars later than their more massive coun- 
terparts. Hierarchical models for galaxy formation produce 
a radically different picture: low mass objects are generally 
predicted to form earlier than more massive ones and con- 
sequently produce older and "redder" stellar populations, 
contrary to observations. In Somerville and Primack (1998), 
dust extinction, assumed to be more important in more mas- 
sive galaxies, is invoked in order to reverse the trend and 
bring agreement between models and observations. In our 
models dust extinction plays a minor role, merely enhanc- 
ing an already correct trend, while at the same time several 
other observables are naturally reproduced. 

The plan of the paper is as follows: In Sec. 2 we present 
the basic ingredients and the underlying assumptions of the 
model. The main input parameters and an application to the 
Milky Way (already discussed in detail in Paper I) are pre- 
sented in Sec. 2.1. The extension to other galactic disks with 
the help of simple scaling relations (established in the frame- 
work of CDM models) is presented in Sec. 2.2. The assump- 
tions about the infall and star formation rate are discussed 
in more details in Sec. 2.3. Some results of the models (evo- 
lution of the gaseous and stellar content, star formation rate, 
colour profiles etc.), as well as the main present day prop- 
erties, are presented in Sec. 3. The results are compared to 
observations in Sec. 4: disk sizes and central surface brigth- 
ness (Sec. 4.1), TuUy-Fisher relations in various wavelength 
bands (Sec. 4.2), colour-colour and colour- magnitude rela- 
tions (Sec. 4.3), gas fractions vs. magnitudes and colours 



(Sec. 4.4), abundances vs. local and integrated properties 
(Sec. 4.5), integrated spectra for different galactic rotational 
velocities (Sec. 4.6). The results are summarised in Sec. 5, 
where some comments are also made on the emerging overall 
picture. 



2 THE MODEL 

In the previous paper of this series (BP99) we presented 
a detailed model of the chemical and spectrophotometric 
evolution of a spiral galaxy and we applied it to the case 
of the Milky Way disk. In this section, we recall briefly the 
main ingredients of the model and the results for our Galaxy 
(Sec. 2.1). Assuming that the Milky Way is a typical spiral 
galaxy, we extend the model to the study of other spirals, 
helped by some simple "scaling laws" (Sec. 2.2), which have 
been established in the framework of the Cold Dark Matter 
scenarios for galaxy formation. The SFR is then calculated 
in a self-consistent way, while a further "scaling law" has 
to be assumed for the infall rate, a crucial ingredient in all 
realistic models of spiral evolution (see Sec. 2.3). 

2.1 Main ingredients and application to the Milky 
Way 

The galactic disk is simulated as an ensemble of concentric, 
independently evolving rings, gradually built up by infall 
of primordial composition. The chemical evolution of each 
zone is followed by solving the appropriate set of integro- 
differential equations, without the Instantaneous Recycling 
Approximation. Stellar yields are from Woosley and Weaver 
(1995) for massive stars and Renzini and Voli (1981) for in- 
termediate mass stars. Fe producing SNIa are included, their 
rate being calculated with the prescription of Matteucci and 
Greggio (1986). The adopted stellar Initial Mass Function 
(IMF) is a multi-slope power-law between 0.1 Mq and 100 
M0 from the work of Kroupa et al. (1993). 

The star formation rate (SFR) is locally given by a 
Schmidt-type law, i.e. proportional to some power of the gas 
surface density Eg: 'I' oc Eg'^, according to the observations 
of Kennicutt (1998); it varies with galactocentric radius R, 
as: 

*(t,i?) = a Eg(t,i?)'-^ V(i?) i?"^ (1) 

where V{R) is the circular velocity at radius R. This radial 
dependence of the SFR is suggested by the theory of star 
formation induced by density waves in spiral galaxies (e.g. 
Wyse and Silk 1989). Since V{R) ~const. in the largest part 
of the disk, this is equivalent to '^{R) oc Eg(ii)^'^ii~^ This 
latter form of the SFR is used in Prantzos and Silk (1998) 
and BP99. However, when extending the model to other 
galaxies with different rotational velocities, we use the form 
of Eq. (1), and we drop the assumption of V{R) — const. 
by calculating the full rotation curve (see Sec. 2.3). 

The infall rate is assumed to be exponentially decreas- 
ing in time, i.e. 

f{t,R) = A{R) e"'/"'-^) (2) 

with t{Ro=8 kpc) = 7 Gyr, in order to reproduce the local 
G-dwarf metallicity distribution and t{R) increasing out- 
wards, from t{R=2 kpc) = l Gyr to t(R=17 kpc)=12 Gyr 
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(see Fig. 3). This radial dependence of f{R) is simulating the 
inside-out formation of galactic disks and, combined with 
the adopted ^'(-R) allows to reproduce the observed current 
profiles of gas, oxygen abundance and SFR in the Milky 
Way disk (see BP99). The coefficient A{R) is obtained by 
the requirement that at time T=13.5 Gyr, the current mass 
profile of the disk E(i?) is obtained, i.e. 



f{t,R) = E(i?) 



(3) 



with E(i?) oc e ^/^^ and a scalelength Rg=2.Q kpc for the 
Milky Way disk. 

The spectrophotometric evolution is followed in a self- 
consistent way, i.e. with the SFR *(t) and metallicity Z{t) 
of every zone determined by the chemical evolution, and 
the same IMF. The stellar lifetimes, evolutionary tracks and 
spectra are metallicity dependent; the first two are from the 
Geneva library (Schaller et al. 1992, Charbonnel et al. 1996) 
and the latter from Lejeune et al (1997). Dust absorption 
is i nclude d according to the prescriptions of Guiderdoni et 
al. (1998) and assuming a "sandwich" configuration for the 
stars and dust layers (Calzetti et al. 1994). 

This model has several ingredients (stellar yields, 
tracks, lifetimes, spectra) calculated by the (presumably well 
understood) theories of stellar atmospheres, evolution and 
nucleosynthesis, while it has relatively few free parameters. 
Indeed, we assume that the IMF (obtained by local observa- 
tions with particular care to the evaluation of the low mass 
part by Kroupa et al. 1993) and the dependence of the SFR 
on gas surface density (derived from observations of exter- 
nal spirals, Kennicutt 1998) are not free parameters. On the 
other hand, the efficiency a of the SFR (Eq. 1) is fixed by the 
requirement that the local gas fraction ag{Ro—8 kpc)~0.2, 
is reproduced at T=13.5 Gyr. We consider then that the 
really "free" parameters of the model are the radial depen- 
dence of the infall timescale r{R) and of the SFR; in the 
latter case, the radial dependence (the V{R)/R factor) has 
some theoretical motivations (see Wyse and Silk 1989). 

The number of observables explained by the model is 
much larger than the number of free parameters. Leav- 
ing aside local observables (age-metallicity relationship, G- 
dwarf metallicity distribution, 0/Fe evolution) the model 
reproduces present day "global" ones (profiles of gas, O/H, 
SFR, as well as supernova rates). Moreover, the current disk 
luminosities in various wavelength bands are successfully ob- 
tained along with the corresponding radial profiles; in par- 
ticular, the adopted inside-out star forming scheme leads to 
a scalelength of ~4 kpc in the B-band and ~2.6 kpc in the 
K-band, in agreement with observations (see BP99 for refer- 
ences and details). It is the first time, to our knowledge, that 
a consistent model for the chemical and photometric evolu- 
tion of the Milky Way is developed. Its success encourages 
us to extend it to other disk galaxies, in the way described 
below. 



2.2 Scaling properties of disk galaxies 

For a simplified description of disk galaxies we adopt here 
the "scaling properties" derived by Mo, Mao and White 
(1998, hereafter MMW98) in the framework of the Cold 
Dark Matter (CDM) scenario for galaxy formation. Accord- 



ing to this scenario, primordial density fiuctuations give rise 
to haloes of non-baryonic dark matter of mass M, within 
which baryonic gas condenses later and forms disks of max- 
imum circular velicity Vc- The mass of the haloes is 



M 



10 G H{z) 



(4) 



where H{z) is the Hubble parameter at the redshift z of halo 
formation and G the gravitational constant. The mass of the 
disk Aid is a fraction md of the halo mass 



Md = rud M 



(5) 



It is assumed that the disk is thin, rotationally supported 
and has an exponential surface density profile 



E(i?) 



(6) 



with central density Eq, scalelength Rd and corresponding 
mass 



Md 



2 TT Eo i?d 



(7) 



Eq. (6) describes the initial gaseous profile of the disk; 
in the terms of our infall model of Sec. 2.1, E(7?) corresponds 
to the right hand side of Eq. (3). Eqs. (4), (5) and (6) allow 
to relate the parameters Eq and Rd to an observable, the 
circular velocity of the disk Vc, but not in a unique way 
(disks of the same mass Md and circular velocity Vc may 
have different Eq and -Rd). In order to break the degeneracy 
one may adopt the spin parameter A as a second parameter, 
as e.g. in MMW98. This parameter is related to the halo 
mass M and angular momentum J through 

A = J G"' M~^/^ (8) 

where the total energy E of the halo (assumed to be an 
isothermal sphere) is given by 



E 



(9) 



Finally, assuming that the angular momentum of the disk 



Jd = 2 Md 14 Rd 

is a fraction jd of that of the halo, i.e. 



■- 3d J 



(10) 



(11) 



one may express the disk parameters Eo and Rd in terms of 
Vc and A, as well as jd, rud and H{z): 



Rd = 
and 
Eo = 



-J— A Vc ) H{z) 

10 V2 \nidJ 



10 ,-2,, ( jd\ 

— — md A Vc — 

TT G \md/ 



H{z) 



(12) 



(13) 



In this work we ignore the effects of the formation time 
of the disk, entering Eqs. (12) and (13) through the term 
H(z). Indeed, it is difficult to incorporate it in the framework 
of our infall model since the disk is slowly built up and the 
term "time of formation" has not a precise meaning. We 
will assume then that all disk galaxies start forming their 
stars at the same time (or redshift) and have today the same 
age as the Milky Way, i.e. T—13.5 Gyr (notice that the 
precise formation time is not very important, as far as the 
corresponding redshift is 2 > 3). We will discuss, however. 



© 1999 RAS, MNRAS 000, 



4 iS". Boissier and N. Prantzos 




the Milky Way's dark halo. Using equs. 5, 8, 9, 10 and 11 
one can also evaluate the spin parameter A as: 



100 200 300 

Velocity (km/s) 

Figure 1. Main properties of our model disks. From top to 
bottom: mass, central surface density and scalelength, respec- 
tively. They are plotted as a function of circular velocity Vq and 
parametrised with the spin parameter A. Disk mass depends only 
on Vc- Filled symbols correspond to the Milky Way model, used 
for the scaling of all other disks (through eq. 14 and 15, with 
md = mac=0.05). 



— I i?d M- 

3d 



(16) 



Using the above numerical values for the Milky Way (and 
m.d/ jd=i), we find Ag=0.03. Eqs. 14 and 15 allow to de- 
scribe the mass profile of a galactic disk in terms of the one 
of our Galaxy and of two parameters: Vc and A. The range 
of observed values for the former parameter is 80-360 km/s, 
where for the latter numerical simulations give values in the 
0.01-0.1 range, the distribution peaking around A ~ 0.04 
(MMW98). Although it is not clear yet whether Vc and A 
are independent quantities, we treat them here as such and 
construct a grid of 25 models caracterised by Vc = 80, 150, 
220, 290, 360 km/s and A = 0.01, 0.03, 0.05, 0.07, 0.09, re- 
spectively. The corresponding main properties of the disks 
(mass Md, central surface density Eq and disk scalelength 
Rd) appear in Fig. 1, expressed in terms of the adopted pa- 
rameters Vc and A. Larger values of Vc correspond to more 
massive disks and larger values of A to more extended ones. 

As we shall see in Sec. 4.1 the resulting disk radii and 
central surface brightness are in excellent agreement with 
observations, except for the case of A=0.01. We could have 
excluded that value from our grid on the basis of instability 
arguments. Indeed, disks dominated by their self-gravity are 
likely to be unstable to the formation of a bar. According to 
the discussion in MMW98, for md=0.05, this would happen 
for disks with A <0.03 (see their Fig. 3). We prefer, however, 
to keep this A value in our grid, for illustration purposes; as 
we shall see, this "unphysical" value leads to the production 
of galaxies ressembling to bulges or ellipticals, rather than 
disks. 



later the possibility that present day disks have different 
ages, which is equivalent to assuming different redshifts for 
their formation. 

It is usually assumed that nid and jd are constants and, 
moreover, that jd/md=l, i.e. that the specific angular mo- 
mentum of the material that forms the disk and the dark 
halo are the same (see MMW98 and references therein). 
The latter assumption, although not supported by numeri- 
cal simulations, simplifies considerably our modelling (other- 
wise, radial inflows should be induced by angular momentum 
transfer in the disk, invalidating our model of independently 
evolving rings). In that case, the profile of a given disk can be 
expressed in terms of the one of our Galaxy (the parameters 
of which are designated hereafter by index G): 



^ A iA 

RdG Ag Vcg 

and 

So ^ _md_ (2lY^ Vc 
EoG mdc VAg/ Vcg 



(14) 



(15) 



where we have: RdG='2-Q kpc, Eog = 1150 M0 pc~ , 
Vgg=220 km/s (see BP99 for references on the properties of 
the Milky Way). Notice that the assumption that all galaxies 
started forming at the same redshift eliminated H{z) from 
Eq. 14 and 15. The mass of the disk is MdG=5 10^" Mq and, 
assuming md=0.05 we obtain Mg = 10^^ M© for the mass of 



2.3 SFR and Infall 

The two main ingredients of the model, namely the Star 
Formation Rate ^{R) and the infall time-scale r{R), are 
affected by the adopted scaling of disk properties in the fol- 
lowing way. 

For the SFR we adopt the prescription of Eq. (1), with 
the same efficiency a as in the case of the Milky Way. In 
order to have an accurate evaluation of V{R) across the 
disk, we calculate it as the sum of the contributions of the 
disk (with the surface density profile of Eq. 6) and of the 
dark halo, with a volume density profile given by: 



p{R) 



Po 



(17) 



where po is the central density and Rc the core radius (pro- 
file on a non-singular isothermal sphere). Other halo profiles 
(e.g. Navarro et al. 1997) are, probably, more realistic, but 
their effect on our results is negligible (the resulting rota- 
tion curves differ by less than 20% from the one adopted 
here except in the innermost regions, where spiral waves are 
less effective and the presence of bulges and bars makes the 
SFR very uncertain, anyway). The corresponding contribu- 
tions of the disk and the halo to the rotational velocity are 
given by (Navarro 1998): 

Vd{R) = 2n GT,o RdX [Io{x)KO(x) - Il{x)Kl{x)]^ (18) 
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Figure 2. Rotation curves of our disk models as funetion of ra- 
dius, expressed in kpc {upper pannel) and in disk scalelengths 
{lower pannel). For each of the five values of our grid Vq (80, 150, 
220, 290 and 360 km/s) we plot the results for the corersponding 
five A values {dotted for 0.01, solid for 0.03, short-dashed for 0.05, 
long-dashed for 0.07 and dot-dashed for 0.09). 



where x — R/Rd and In and Kn are modified Bessel func- 
tions of order n ; and 



o Rn R 
Vh{R) = ^Tv G po rI[1- ^arctan{ — ] 

K Rc 



(19) 



respectively, while the final rotation curve is calculated by: 



ViRf 



VhiRf + V4RY 



(20) 



The calculated rotation curves for our models appear in Fig. 
2, as a function of radius expressed in kpc {upper pannel) 
and in scalelengths Rd {lower pannel); it can be seen that 
in the latter kind of "universal" rotation curve is 

obtained (excpet for the unrealistic, as we shall see below, 
A=0.01 case). Although this velocity profile is not exactly 
the same as the one derived observationally in Salucci and 
Persic (1997), it is very close to it and certainly sufficient 
for the purposes of this work. It allows to calculate the SFR 
(Eq. 1) in a fully self-consistent way. 

The choice of the infall rate t{R) is more problematic 
than the one of the SFR because there are no theoretical 



principles to guide us. MoUa et al. (1999) adopt an expo- 
nential radial dependence (with a scale-length of 4 kpc) 
for t{R) in the Milky Way disk and a characteristic time 
To = TQ{Md/MdG)~^^^ for the other galactic disks; this scal- 
ing is made on the basis of the association R^ff/Ro (effec- 
tive radius/solar galactocentric distance) between the Milky 
Way and the galaxies of their study. 

In our case, we have no a priori idea on which region 
of a given disk will be equivalent to the solar neighborhood. 
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Figure 3. Adopted relation between infall time scale r and total 
surface density S, generalized from the Milky Way model {solid 
curve) to other galaxies (see Sec. 2.3). Infall is more rapid in 
denser regions and in more massive galaxies. Negative values of 
l/r mean that the infall rate is increasing in time. 



We then adopt the following prescription: the infall time 
scale increases with both surface density (i.e. the denser in- 
ner zones are formed more rapidly) and with galaxy's mass, 
i.e. T[Md,'S{R)]. In both cases it is the larger gravitational 
potential that induces a more rapid infall. Our adopted infall 
rate is given by: 



f[Md,J:{R),t] oc exp 



t 



^g(E) 



+ 0.4t 1 



Vc{Md 
220 



(21) 



The corresponding characteristic infall timescales 
T[Md, S(ii)] appear in Fig. 3, as a function of E(ii) and 
of Vc- The radial dependence of r on S(i?) is calibrated 
on the Milky Way (figure |^ where r is fixed to 7 Gyr at 
Ro~8 kpc (to solve the G-dwarf problem) and to 1 Gyr at 
2kpc from the galactic center. On the other hand, the mass 
dependence of r is adjusted as to reproduce several of the 
properties of the galactic disks that will be discussed in Sec. 
4. The adopted prescription allows to keep some simple scal- 
ing relations for the infall in our models and the number of 
free parameters as small as possible. We recognise that a 
different parametrisation could be made with a - certainly - 
non negligible impact on the results. In fact, by suppressing 
the dependence of r on galaxy mass (i.e. by assuming the 
same r[E(_R)] relation for all galaxies) we found that several 
properties of disks - in particular, the colour magnitude and 
metallicity-luminosity relations - are not well reproduced. 
The fact that the adopted simple prescription provides a 
satisfactory agreement with several observed relationships 
in spirals (see Sec. 4) makes us feel that it could be ulti- 
mately justified by theory or numerical simulations of disk 
formation. 



3 MODEL RESULTS 

We have computed the chemical and specrophotometric evo- 
lution of a grid of 25 disk models (defined by the values of 
Vc and A given in Sec. 2.2) up to an age of T=13.5 Gyr. We 
calculated the evolution of the various chemical and pho- 
tometric profiles, as well as the disk integrated properties 
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Figure 4. Histories of total SFR solide curve), gas mass (g, short dashed curve) and stellar mass (s, long dashed curve) for our 
disk models. Each column corresponds to a value of A and each row to a value of Vq- The dotted curve is an exponential fit to the SFR. 



according to the procedure outlined in paper I (BP99, Sec. 
3.3). 

In this work, we concentrate mainly on the integrated 
properties at T=13.5 Gyr and compare them to observations 
of local spirals (i.e. at redshift z ~0, see Sec. 4). However, in 
this section, we also present some of the results concerning 
photometric and colour profiles, since their analysis helps to 
understand the behaviour of the integrated properties and 
the role of the adopted prescriptions for the SFR and the 
infall rate. The implications of our models for galaxy obser- 
vations at higher redshifts, and the corresponding chemical 
and photometric gradients will be analysed in forthcoming 
papers. 



3.1 Histories of gas, stars and SFR 

The evolution of the SFR and of the gaseous and stel- 
lar masses in our 25 disk models is shown in Fig. 4. This 
evolution results from the prescriptions for infall and SFR 
adopted in Sec. 2.3. Values of SFR range from ~100 Mq/jt 
in the early phases of massive disks to ~0.1 Mq/jt for 
the lowest mass galaxies. The resulting SFR history is par- 
ticularly interesting when compared to simple models of 
photometric evolution, that are usually applied in studies 
of galaxy evolution and of their cosmological implications. 
Such models (one-zone, no chemical evolution considered in 
general) adopt exponentially declining SFR with different 
caracteristic timescales for each galaxy morphological type 
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Figure 5. Caracteristic time-scales for star formation, i.e. for 
forming the first half of the stars of a given galaxy, resulting from 
our models. They are plotted as a function of the spin parameter A 
[left) and of circular velocity Vc (right). Filled points correspond 
to the Milky Way. 



(Guiderdoni & Rocca-Volmerang e 198 7, Bruzual & Chariot 
1993, Fioc & Rocca-Volmerange ([l997| )). 

An inspection of Fig. 4 shows that this may be a rela- 
tively good approximation for a large part of galaxy's his- 
tory, but not for the earliest and latest times (i.e. for the first 
and last couple of Gyr) . We used an exponentially declining 
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Figure 6. Evolution of the B-band surface brightness profile of our disk models. Profiles are shown at ages T=3 (dotted), 7.5 (dashed) 
and 13.5 Gyr (solid). The grey horizontal line at fig=22.5 mag arcsec"^ indicates the limit for Low Surface Brightness galaxies. Each 
column corresponds to a value of A and each row to a value of Vc (in km/s). 



SFR to fit the obtained SFR histories (shown by the dot- 
ted curves in Fig. 4). The shortest caracteristic timescales 
correspond to the more massive and compact disks (larger 
Vc, smaller A): in those cases, the rapid infall makes the 
gas available early on, while the large Vc/Rd factor of the 
adopted SFR leads to the rapid consumption of that gas. 
Whether such massive and compact disks exist is another 
matter (see Sec. 4). As we move to larger values of A (i.e. 
more extended disks with larger Rd), the Vc/Rd factor 
makes star formation less efficient and the resulting SFR has 
a smoother history. In fact, for most of the values of the pa- 
rameter couple (Vc, A) the SFR remains essentially constant 
during most of the disk's history, since the gas consumption 
is compensated by the infall. Finally, in the case of the most 
extended and less massive galaxies (low Vc, large A) the SFR 
and gas mass actually increase in time, due to the adopted 
prescription for infall (which is very slow for low mass and 
low surface density galaxies). Extended galaxies are found to 
have low surface brightness (as already shown in Jimenez at 
al. 1998) and high gas fractions (see Sec. 3.2). In summary, 
it turns out that A affects mainly the shape of star formation 
history (i.e. the caracteristic time scale of star formation), 
while Vc determines the absolute values of SFR, gaseous and 
stellar content (see Sec. 3.4). 

The fact that in some cases the SFR increases in time 
leads us to define the time-scale r, , required for forming the 
first half of the stars of a given galaxy. This time scale for 
star formation appears in Fig. 5, plotted as a function of A 



(left) and Vc (right). It can be seen that r, is a monoton- 
ically increasing function of A and a decreasing function of 
Vc- For the Milky Way (Vcg=220 km/s, Ag=0.03) we find 
r» ~3.5 Gyr. For galaxies less massive than the Milky Way 
and A > 0.06, it takes more than half the age of the Universe 
to form the first half of their stars. 

3.2 Photometric profiles 

In this section we present some results on the photometric 
evolution of our disk models. We illustrate this evolution by 
commenting only on the B-band surface brightness profile. 
Other photometric bands present similar trends and are not 
shown here, for lack of space; results concerning some of 
them will be presented in Sec. 4, where comparison is made 
to a large body of observational data. 

The evolution of the B-band surface brightness is shown 
in Fig. 6 for three ages: t=3, 7.5 and 13.5 Gyr, respectively. It 
reflects clearly the adopted scheme of inside-out star forma- 
tion, i.e. the fact that gas is more rapidly available in the in- 
ner galactic regions (because of the adopted infall timescale 
t(R)), where it is efficiently turned into stars (because of 
the V/R dependence of the SFR). At late times, relatively 
little gas is left in the inner regions, but quite a lot is still 
infalling in the outer ones, maintaining an important SFR 
there. Several interesting features should be noticed: 

a) The central surface brightness hbo today depends 
little on Vc, while the more compact (higher A) disks have. 
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Figure 7. Evolution of the B-K profile of our disk models. Profiles are shown at ages T=3 (dotted), 7.5 (dashed) and 13.5 Gyr (solid). 
Each column corresponds to a value of A and each row to a value of Vq (in km/s). Colour profiles generally fiatten with time, because 
of the inside-out scheme of star formation adopted (see text). 



in general, higher /iso. In almost all cases fiBO was higher 
in the past 

b) As a result of the inside-out star formation, the 
B-band surface brightness profile flattens, in general, with 
time. Disks were more compact in the past. 

c) Photometric profiles are rarely "purely exponential" ; 
indeed, they may deviate considerably from the usually 
adopted exponential ones. This is in agreement with the 
conclusions of Courteau and Rix (1999) who find that only 
about 20% of their sample of more than 500 disk galaxies 
qualify as pure exponentials. 

d) Because of the inside-out star formation scheme, 
the isophotal radius R25 (the radius at which /is =25 
mag arsec"^) generally increases with time (except for the 
A=0.01 case). Present day values of R25 in our models are 
found to range from a few kpc to ~40 kpc. 

e) Several of the models in Fig. 6 have hbo <22.5 
mag arcsec"^ at present, i.e.they correspond to Low Surface 
Brightness (LSB) galaxies. This concerns all the galaxies 
with A=0.09 and some with A=0.07. However, according to 
our models, some of those galaxies have always been LSBs 
and some of them not. In the former class belong the less 
massive galaxies, which were fainter in the past because of 
their continuously increasing SFR (Fig. 4) . In the latter class 
belong more massive galaxies, which have a declining SFR 
and therefore were brighter in the past; they became LSBs 
only after a certain age. Finally, in some intermediate cases 
(e.g. Vb=150 km/s and A=0.05, 0.07) the central SFR and 



corresponding fiBo have remained virtually unchanged dur- 
ing the galaxy's lifetime. These results may have important 
implications for the study of LSBs as we shall see in Sec. 
3.4. 

3.3 Colour profiles 

Figure 7 presents the evolution of the B-K profiles of our 
model disks. Again, the SFR histories of Fig. 4 and the 
inside-our star formation scheme allow to understand the 
main features: 

i) The more compact and more massive a galaxy is, the 
larger is the final B-K value, either the central (Fig. 7) or 
the integrated one (Fig. 8). This is due to the fact that low 
A and high Vc values favour a rapid early star formation 
and the existence of relatively old stellar populations today. 
In fact, the properties of A=0.01 galaxies correspond more 
to those of bulges or ellipticals than disk galaxies. 

ii) B-K values are always larger in the inner galactic re- 
gions, since stars are formed earlier there than in the outer 
disks. As time goes on, a "reddening" wave propagates out- 
wards in the disks. 

iii) More massive disks have a flatter B-K profile than 
less massive ones with the same A. 

iv) The outer galactic regions have, in general, a "flat" 
B-K profile, because they form stars at late times. Stellar 
populations in those regions have, on the average, the same 
age (they are young) and same B-K values. On the other 



© 1999 RAS, MNRAS 000, 



Chemo-spectrophotometric evolution of spiral galaxies:!!. Main properties of present day disk galaxies 



hand, in the inner regions there are substantial discrepancies 
in the timescales of star formation, and colour gradients are 
consequently established. 

v) For a given Vc the final B-K profile is flatter in more 
compact disks, since the 1/R factor of the SFR and the infall 
timescales are such that efficient early star formation may 
take place even in their outer regions. 

The results presented in the previous section and in 
this one may have interesting implications for LSB galax- 
ies. Indeed, it can be seen in Figs. 5 and 6 that disks 
with the same fiBO today may have quite different central 
colours. This is e.g. the case for (A = 0.07, Vc=80 km/s) and 
(A = 0.07, Vc=220 km/s). They both have hbo ~23 mag 
arscec"^, but the former has (B-K)o=2.8, while the latter 
has (B-K)o=4. This is due to their different SFR histories, 
since in the former case most of the stars are formed late, 
while in the latter they are formed quite early on. These re- 
sults could be related to the conclusions of a recent study by 
Bell et al. (1999), who find that there are two kinds of LSBs, 
"red" and "blue" , with different star formation histories: the 
former are better described by a "faded disk" senario, while 
the latter are described by models with low, roughly con- 
stant star formation rate. According to our senario, these 
different routes are dictated simply by the galaxy's mass, 
which induces an earlier SFR in the former case (because of 
the more rapid infall). A clear prediction of our models is 
that "blue" LSBs should be less massive and less luminous 
than the "red" ones. 

The observational implications of colour and photom- 
etry profiles and their relationship to abundance gradients 
(not presented here) is postponed to a forthcoming paper 
(Boissier and Prantzos, in preparation). 

3.4 Integrated present day properties 

Assuming that all our model disks started forming their stars 
~13 Gyr ago, we present in Fig. 8 their integrated properties 
today (at redshift z—0). Several interesting results can be 
pointed out: 

1) For a given A, the stellar and gaseous masses are 
monotonically increasing functions of Vc- For a given Vc 
disks may have considerably different gaseous and stellar 
contents (by factors 30 and 5, respectively), depending on 
how "compact" they are, i.e. on their A value. The more 
compact galaxies have always a smaller gaseous content. 
This is well illustrated by the behaviour of the gas fraction 
a g {bottom left of Fig. 8). There seems to be a rather good 
correlation between ag and Vc for all A values, a feature that 
is indeed observed (see Sec. 4.4). The most compact galaxies 
have always ag lower than 0.1, while the most extended ones 
(those corresponding to the LSBs discussed in the previous 
section) have always ag >0.5. The latter property seems in- 
deed to be an observed feature of LSBs (McGaugh and de 
Blok 1997, and Sec. 4.4). In general, the final gas fraction is 
a slowly decreasing function of Vc and a rapidly increasing 
function of A. 

2) The pattern of the average oxygen abundance (i.e. 
total oxygen mass in the gas, divided by the gas mass, central 
pannel in Fig. 8) can be readily understood on the basis of 
the behaviour of the gas fraction: wherever the latter is low 
(see point 1) the metallicity is high. Most of our models have 
average metallicities lower than solar. Only the most massive 



and compact ones develop average metallicities larger than 
2 times solar. Our LSB candidates have average metallicities 
~0.5 solar. 

3) The total SFR presents an interesting behaviour 
{middle bottom pannel): the current SFR is found to 
be mostly a function of Vc alone, despite the fact that the 
amount of gas varies by a factor of ~2-3 for a given Vc- 
This is due to the fact that more extended disks (higher A) 
are less efficient in forming stars, because of the 1/Rd fac- 
tor (for a given Vc)- This smaller efficiency compensates for 
the larger gas mass available in those galaxies, and a unique 
value is obtained for 'i'o- This immediately implies that little 
scatter is to be expected in the B-magnitude vs. Vc relation, 
and this is indeed the case as can be seen in Fig. 8 {middle 
right pannel)- The A=0.01 case does not follow this pattern, 
since very little gas is left in the disk; as stressed on several 
occasions, this case corresponds rather to bulges or ellipti- 
cal galaxies and is shown here only as a limiting case of our 
disk models. The implications for the TuUy-Fisher relation 
are discussed below (Sec. 4.2). 

4) The mass to light ratio T = AI/L is a quantity 
often used to probe the age and/or the IMF of the stel- 
lar population and the amount of dark matter, but also 
to convert results of dynamical models (i.e. mass) into ob- 
served quantities (i.e. light in various wavelength bands). 
From disk dynamics of high surface brightness galaxies Bot- 
tema (1997) inferred Ts=1.79±0.48 (for a Hubble parame- 
ter /i=0.75). Using the value of B — 1=1.7 for typical disk 
galaxies given by McGaugh and de Blok (1997), MMW98 
derived T/ = (1.7±0.5)/i, where the Hubble parameter Ho = 
100 h km/s/Mpc. 

Our model values for Tb and T/ are plotted in Fig. 8 
{middle upper pannel) along with the corresponding obser- 
vational range {grey bands)- It can be seen that: a) both Tb 
and Ti model values are slowly decreasing functions of Vc 
and slowly increasing functions of A; b) most of the model 
values lie well inside the observationnaly derived range; c) 
the A=0.01 case lies outside the observed range of Tg (an- 
other indication that this value of the spin parameter does 
not produce realistic disks). 

In summary, our models lead to fairly acceptable values 
of Ts and T/, these values are confined to a relatively nar- 
row range and a systematic (albeit weak) trend is obtained 
with galaxy's mass, that should not be neglected in detailed 
models. 

5) As discussed in Sec. 3.3, the central surface bright- 
ness fiBO depends little on Vc, but quite strongly on A. This 
can be qualitatively understood from Eq. (15) implying for 
the central surface density a strong dependence on A, but a 
mild one on Vc- Taking into account that the A-distribution 
peaks around 0.04-0.05 (MMW98) we see that most of our 
model disks have a roughly constant current jiBO ~21.5 mag 
arcsec"'^. This is quite an encouraging result, in view of the 
observed constancy of the central surface brightness in disks, 
around the "Freeman value" /iso=21.7±0.3 mag arcsec"^ 
(Freeman 1970); this value is indicated by the grey band in 
Fig. 8 {top right pannel)- In the next section we show that 
the obtained values of central surface brightness in longer 
wavelengths are also in fair agreement with recent observa- 
tional results. Finally, models with A=0.07 and 0.09 result 
in LSB galaxies today (but not always in the past, see the 
discussion in Sec. 3.3). 
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Figure 8. Integrated properties of our models at an age T=13.5 Gyr. Mass of stars {top left), mass of gas (middle left), gas fraction 
{bottom left), M/L values in B- and I- bands {top center), average oxygen abundance, i.e. total mass of oxygen divided by total mass 
of gas in the disk {middle center), star formation rate {bottom center), central surface brightness fiBO {top right), total B-magnitude 
(middle right), total B-K (bottom right). All results are plotted as a function of the rotational velocity Vc and parametrised with the 
values of the spin parameter A. The filled symbols on the A=0.03 curves and for Vc'=220 km/s in all pannels correspond to Milky Way 
values; as stressed in Sec. 2 the Milky Way is used to calibrate all our models, through the adopted scaling relations (eqs. 14 and 15). 
The grey bands in the M/L and fiBO pannels indicate the corresponding ranges of observed values (see Sec. 3.4, points 4 and 5). 



6) The final B-K colours of our models (Fig. 8, bottom 
right pannel) can also be readily understood by an inspection 
of the gas fraction pattern: for a given A, more massive disks 
have less gas and are redder, while for a given Vc more ex- 
tended disks (larger A) have more gas and are bluer. Overall, 
the less massive galaxies have B-K values in the 2.5-3 range, 
while for the most massive ones we obtain values around 
3.5-4. 



4 RESULTS VS OBSERVATIONS 

Having described in the previous section the main results of 
our models in terms of the input parameters Vc and A, we 
turn now to a comparison of these results to observations. 
Our purpose is to check whether the results of our grid of 
models (a) bracket reasonably well the range of observed 
values of various galactic properties and (b) reproduce some 
observed strong correlations. Such a comparison is, in prin- 
ciple, possible if the age of the observed galaxies is known, 
but this is unfortunately not the case. We shall make then 
the assumption that all galaxies started forming their stars 
at the same epoch, around 13 Gyr ago. It is not a priori 
obvious whether this assumption is reasonable or not, but 
we shall see that in most cases it leads to quite acceptable 



results. In any case, as stressed in Sec. 3.1 the term "time 
of formation" is a rather loose concept for models with in- 
fall: the "effective" timescales for the formation of the stellar 
populations in our models are given in Fig. 5. 



4.1 Boundary conditions: Disk size and central 
surface brightness 

As explained in Sec. 2.2, our model disks are described by 
the central surface density Eo and scalelength Rd- These are 
related to the "fundamental" parameters Vc and A (fun- 
damental because they refer to the dark matter haloes) 
through the scaling relations (14) and (15), which allow to 
use the Milky Way as a "calibrator" . 

In this section we check whether this type of modelling 
produces present day disks with "reasonable" morphological 
properties. The central surface brightness and scalelength in 
the r-band represent fairly well the underlying stellar pop- 
ulation and are less affected by uncertainties due to dust 
extinction. The behaviour of our models in shorter wave- 
lengths is somehat different, but we postpone an analysis of 
the colour gradients of our models to a forthcoming paper. 

In a recent paper, Courteau and Rix (1999) analyse two 
large data samples of more than 500 non-barred disk galax- 
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Giovanelli et al. 
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TuUy et al. 


8.17 


-21.54 


0.42 mag 


Han & Mould 


7.87 


-21.42 


0.40 mag 


K^atliewsoii 


6.80 


-21.72 


0.43 mag 


Our models 








A=0.01 


6.81/7.17 


-21.34/-21.70 




A=0.03 


7.01/7.31 


-21.47/-21.72 




A=0.05 


7.37/7.71 


-21.35/-21.54 




A=0.07 


7.79/8.10 


-21.16/-21.30 




A=0.09 


8.18/8.44 


-21.00/-21.05 





Table 1. The TuUy-Fislier relation in the I-band: observations 
vs. our models (see Sec. 4.2 for references). A value of /i=0.65 
is assumed for the Hubble parameter in this Table and Fig. 10. 
Two values (separated by /) are given for our models: the first one 
corresponds to models with extinction taken into account and the 
second to models with extinction neglected. Observational data 
are from Giovanelli et al. (1997), TuUy et al. (1998), whereas the 
Han and Mould and Mathewson et al. data are taken from the 
compilation of Willick et al. (1996). 



Figure 9. Central surface brightness (upper pannel) and disk 
size [lower pannel) in the r-band, as a function of the galaxy's 
magnitude M^. Our results are parametrised by the values of the 
spin parameter A; A=0.01 models produce unrealistically small 
disks with very high surface brightness and correspond rather to 
galactic bulges. Data are from Courteau and Rix (1999). 



ies. They find that only ~20% of these disks can be described 
by purely exponential profiles and they provide disk sizes, 
rotational velocities, magnitudes and colours. Their data al- 
low to derive the central disk surface brightness. 

In Fig. 9 we present the r-band central surface bright- 
ness (upper pannel) and size (lower pannel) of our model 
disks as a function of the galaxy's luminosity and we com- 
pare them to the data of Courteau and Rix (1999). It can 
be seen that: 

i) The central surface brightness in our models compares 
fairly well to the data. For a given value of A, Uro is quasi- 
independent on luminosity (or Vc), as explained in Sec. 3.4. 
The obtained values depend on A, the more compact disks 
being brighter. All the observed values are lower than our 
results of A=0.01 models. Higher A values bracket reasonably 
well the data. 

ii) The size of our model disks decreases with decreasing 
luminosity, with a slope which matches again the observa- 
tions fairly well. Our disk sizes range from Rd ~6-10 kpc for 
the most luminous galaxies to 1-2 kpc for the less luminous 
ones. The A=0.01 curve clearly lies too low w.r.t. all data 
points, showing once more that it does not produce realistic 
disks. 

It should be noticed that a similar comparison between 
theory and observations is made in MMW98. However, for 
lack of an evolutionary model, MMW98 compare rather the 
profiles of their gaseous disks (i.e. those given by our eq. 12 
and 13) to the data, making the implicit assumption that 
gas is transformed into stars with the same efficiency in all 
parts of all disks. Our models show that this is a good ap- 
proximation for the inner regions, but not for the outer ones 



(because of the adopted radial dependence of the SFR effi- 
ciency) and that the final stellar disk sizes are shorter than 
the initial gaseous ones. Moreover, MMW98 adopt a M/L 
ratio (the same for all disks) in order to convert their mass 
surface density to surface brightness. As discussed in point 
4 of Sec. 3.4, the T = M/L ratio is a slowly increasing func- 
tion of a galaxy's mass (at least accordingly to our models). 
Our self-consistent chemical and spectrophotometric evolu- 
tion models allow to avoid these approximations and to com- 
pare directly the output to observations. 

4.2 Dynamics: the Tully-Fisher relation 

The Tully-Fisher (TF) relation is a strong correlation ex- 
hibited by the whole population of disk galaxies. It relates 
their luminosity L (or magnitude M) and their circular ve- 
locity Vc (or Ha line width W) and can be expressed as 
M= -a [log{2Vc) - 2.5] -I- 6 or L oc Vc^'^''^. The exact slope 
of the relation and its zero point are still the subject of con- 
siderable debate (e.g. Giovanelli et al. 1997). For this reason, 
we compare our results with several recent fits to the I-band 
TuUy-Fischer relation obtained with various data samples. 

In a recent study, Giovanelli et al. (1997) compiled a ho- 
mogeneous set of I-band photometry and HI velocity profiles 
of 555 spiral galaxies in 24 clusters with redshift extending 
to cz ~9000 km/s. The slope of their TF relation is 7.68 ± 

0. 13. They estimate the mean amplitude of the scatter in 
the observed relation to be 0.35 mag, but they notice that it 
increases with decreasing velocity. They also estimate that 
the combination of statistical and systematic uncertainties 
can affect the zero point of the relation up to 0.07 mag, i.e. 
a rather small effect. In order to compare our results to the 
Giovanelli et al (1997) sample we shall assume that h — 0.65 

1. e. a value very close to the h = 0.69 ± 0.05 value derived 
by these authors on th e basis of their data. 

Willick et al.( 1996 ) have computed several TF relations 
on the basis of samples from various authors. Two of those 
samples are in the I-band and shall be used in the present 
work. The first one is a cluster sample from Han, Mould 
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Figure 10. The Tully-Fisher relationship in the I-band. In each 
pannel, our results (curves parametrised by the corresponding A 
value, as indicated in the bottom right pannel) are compared with 
fits from various authors (shaded zones), differing by their zero- 
points and slopes (see Table 1 and Sec. 4.2 for references on the 
data). The models shown here are those with extinction taken 
into account (first column in Table 1), which lead to flatter TP 
relations than when extinction is neglected. In all cases (except for 
TuUy et al. 1998) the data are plotted assuming that the Hubble 
parameter h=0.65. 



and coworkers (references are given in Willick et al.(1996)). 
They find a slope a — 7.87± 0.16 and a scatter of 0.4 mag. 
The second one is the field galaxies sample of Mathewson 
et al. ( 1992| ), composed of 1219 galaxies. It provides a slope 
a = 6.80± 0.08 an d a sca tter of 0.43 mag. For these two sam- 
ples Willick et al.(1996) define the "Hubble absolute mag- 
nitude" of a galaxy, computed from its apparent magnitude 
and redshift and independent of the Hubble parameter; the 
zero-points of the TF relation are then determined through 
this Hubble absolute magnitude. In order to plot those two 
relations on a common scale with other samples, we trans- 
formed the quoted "Hubble absolute magnitudes" into usual 
absolute magnitudes by assumi ng h= 0.65. 

A fourth work (TuUy et al.(1998)) uses a limited sample 



of 87 galaxies from the Ursa Major and Pisces Clusters. The 
obtained slope is steeper than in the other determinations 
(a=8.17) while the scatter is similar (0.42 mag). In that 
work, TF relations in other photometric bands (B,R,I,K') 
are considered for the same galaxy sample; we shall compare 
them to our results in the end of this section. 

The parameters of those various fits are presented in 
Table |l], and the data are compared with our models in Fig. 
10. It can be seen that: 

a) The model slopes and zero-points depend somewhat 



on A: in general, larger A values lead to steeper M-Vc rela- 
tions. This is is due to the following reason: At high Vc all 
galaxies transform most of their gas into stars, i.e. the final 
stellar mass depends little on A and so is the case for the 
corresponding M/. As we go to lower Vc, the final stellar 
mass depends more and more on A, because higher A galax- 
ies are more inefficient in forming stars: this is refiected in 
a steeper slope of the corresponding M-Vc relation. 

b) Since the most probable values of A are in the 0.04- 
0.05 range (according to numerical simulations), our corre- 
sponding slopes of the M-Vc relation are a ~7. 1-7.3 (with 
extinction taken into account). This value of a is intermedi- 
ate between those obtained from the samples of Giovanelli 
et al., Mathewson et al. and Han, Mould and collaborators. 
It is certainly lower than the slope in the TuUy et al. sample, 
which concerns, however, disk galaxies in clusters; in such 
environments disks probably evolve quite differently than in 
the field, being affected by interactions. Our models are bet- 
ter compared to field galaxies, such as those in the sample 
of Mathewson et al. 

c) Because of the very inefficient star formation at low 
Vc and high A, the TF relation of our models shows a mag- 
nitude dependent scatter, from -^0.1 mag at high Vc to 
~0.5 mag at low Vc- This trend is not apparent in the ob- 
served relations appearing in Fig. 10 [shadowed zones), since 
only the average scatter appears. However, Giovanelli et al. 
(1997) find that there is indeed a trend of increasing scatter 
with decreasing magnitude in the TF relation. Notice that 
in our models, the galaxies that contribute mostly to that 
scatter at low Vc are the low surface brightness galaxies 
(A ~0. 07-0. 09); such galaxies are, in general, not included 
in the observed samples. In any case, a clear prediction of 
our models is that the scatter in the TF relation should be 
more important at low Vc- 

d) Extinction modifies somewhat the slopes, leading to 
slightly swallower TF relations than those obtained by the 
stellar population alone. This is because it affects more the 
massive (high surface brightness and high metallicity) disks, 
rather than the low mass ones. If our extinction values are 
overestimated, then the slopes of our models should be in- 
termediate between the two values given in Table 1. This is 
also true for the zero-points. 

It should be noticed here that, if the initial gas were al- 
most completely turned into stars, one should naturally ex- 
pect a. L cc Vc relation (from Eq. 4, assuming — const 
and M/L = const). In other words, the adopted scaling re- 
lations lead to a "built-in" TF relation, as properly pointed 
out in MMW98. However, when star formation is properly 
taken into account in "realistic models" , different slopes are 
found, depending on the SFR efficiency (itself a function of 
A, see Table 1). Although the differences introduced by the 
more realistic senarios are relatively small, they may have 
important implications for a proper interpretation of the TF 
relation. In particular, our models, taken at face value, sug- 
gest that low surface brightnesss galaxies (i.e. those with 
A >0.07) should have steeper TF relations than those of 
high surface brightness in the I-band. 

TuUy et al. (1998) also investigated the slope of the TF 
relation in several wavelength bands. They found that the 
slope increases with wavelength, as can be seen in Table 
2. Our results are also shown on Table 2 and in Fig. 11, 
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Figure 11. The TuUy-Fisher relationship in several photometric 
bands. Our results [curves parametrised by the corresponding A 
value, as in dicate d in the bottom pannel) are compared with the 
Tully et al.(1996) observations (shaded zones). 
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Table 2. Slopes in the TF relations Lx oc V^^'^'^, in various wave- 
length bands. Our models (arranged by A values) are compared 
to the observations of Tully et al. (1998); see also Fig. 11. 



compared to the observational fits of Tully et al. (1998). 
The following points should be noticed: 

i) In our models, the trend of increasing slope with in- 
creasing A is found in all wavelength bands, not only in the 
I- band. 

ii) As we go from the B- to the K'-band, our models do 
predict a steeper and steeper slope for all A values. Leaving 
aside the A=0.01 case (unrealistic for disks) one sees that 
the slope a increases between B- and K'-bands by ~1.8 for 
all A values. This increase is somewhat larger than the cor- 
ersponding increase of ~1 in the data of Tully et al. (1998). 

iii) Our model values for the slope a are, in general, 
smaller than those of Tully et al. (1998), especially in the 
B-band. Only the A=0.09 models lead to values comparable 
to those of Tully et al. (1998) in the R, I and K' bands. 

As stressed already, the Tully et al. (1998) data con- 
cern spirals in clusters which cannot be directly compared to 
those modelled in our work; our slopes in the I-band compare 



rather well to the field galaxy sample of Mathewson et al. 
But the important issue here is the trend of increasing slope 
with wavelength, found in both theory and observations. In 
our models this increase is due to the fact that less massive 
galaxies are chemically and photometrically "younger" than 
their more massive counterparts, having in general "bluer" 
colours. This is shown in Fig. 8 (bottom right pannel) , where 
B-K increases systematically with Vc- This change in B-K 
with Vc produces a steepening in the slope of the TF re- 
lation as we go from the B- to the K-band. This feature 
is not "built-in" in our models (it does not stem directly 
from the adopted scaling relations) but it results from the 
adopted prescription for the SFR and the infall rate (Sec. 
2.3), making less massive galaxies appear "younger" than 
more massive ones. This generic feature of our models is 
crucial also in interpreting the metallicity-luminosity rela- 
tion of disks (see Sec. 4.5) and we consider it to be one of 
the most important points of this work. 

In a recent paper, Heavens and Jimenez (1999) use also 
semi-analytic disk models in CDM senarios and a rather sim- 
ple prescription for star formation in order to investigate the 
variation of the slope of the TF relation with wavelength. 
They apparently manage to reproduce the oberved trend, 
but they offer no explanation for their success. They also 
suggest that A value and formation redshift z are the pa- 
rameters that can explain the observed dispersion in the TF 
relation. Our models show that this is certainly true for A, 
but the disk baryonic fraction nid (taken as constant here) 
could obviously play a similar role. 

4.3 Luminosity evolution: Colour-colour and 
colour-magnitude relations 

Spectro-photometric models are usually calibrated on 
colour-colour and colour-magnitude diagrams. The param- 
eters of the models (star formation efficiency, stellar IMF, 
infall time-scales, ages ...) are chosen as to reproduce the 
propertie s of g alaxies of all morphological t ypes (e.g. Lind- 
ner et al.( 1999| ), Fioc & Rocca- Volmerange( 1997 )) . 

As explained in the Introduction, our study is limited to 
the case of spiral galaxies, because they offer a much larger 
amount of observational data to compare with. However, we 
have not "tuned" any parameters in order to reproduce any 
diagram of this kind. It is thus of interest to see if our mod- 
els, generalized from the Milky Way by the use of simple 
scaling relations, are able to reproduce those diagrams for 
spiral galaxies. For this aim, we used data from De Jong, 
1996. De Jong's data consists in radial and integrated pho- 
tometry of 86 undisturbed face-on spirals. They are thus 
directly comparable to our models, since uncertainties due 
to inclination, or to perturbations by nearby galaxies (which 
could have affected the star formation activity) are avoided. 
We evaluated the disk magnitudes and colours in de Jong's 
sample by using his fits for the disks only (i.e. central sur- 
face brightness and disk scalelength) and plot them on Fig. 
12 (notice that the error bars refiect the uncertainties of the 
fit, not those of the observations). 

Since the age of the observed galaxies is not known, we 
show in Fig. 12 our results for three different epochs: t=3, 
7.5 and 13.5 Gyr, covering the expected range of galaxy ages. 
It is clearly seen that: 

1) At all ages there is a tight B — K vs B — V corre- 
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Figure 12. Colour-magnitude and colour-colour relations for our models at three different ages: 3 Gyr {upper pannels), 7.5 Gyr (middle 
pannels) and 13.5 Gyr {lower pannels). The grid of models represents 5 values for the spin parameter A (0.01 to 0.09 from top to bottom 
in solid curves) and 5 values for the circular velocity Vc (80 to 360 km/s from left to right, in dashed curves), arranged as indicated in 
the lower left pannel. In the B — K vs B — V relation our grid is "squeezed" onto a relatively tight relation, determined mainly by Vq- 
Data in all pannels (the same points are plotted in each row) are from de Jong (1996). We have calculated magnitudes and colours for 
disks only, using de Jong's fits (central surface brightness and disk scalelength). Clearly, the observed galaxies are older than 3 Gyr, and 
rather closer to 13.5 than to 7.5 Gyr. 



lation, quasi-independent of the value of A. It depends only 
on the galaxy's mass (i.e. luminosity): lower mass galaxies 
are always "bluer" than more massive ones. 

2) Results at t=3 Gyr do not compare well to the data: 
model disks are too blue. The observed galaxies cannot be 
as be as young as 3 Gyr. 

3) At t=7.5 Gyr the results compare much better to the 
data, although there is still a lack of "red" model galaxies 
(i.e. with B — V in the 0.8-1.1 range). 

4) At t=13.5 Gyr the agreement between theory and 
observations is quite satisfactory. The observed colours and 
magnitudes are well bracketed by our grid of results. 

5) Our models do predict a colour-magnitude relation 
which becomes steeper with time. More luminous galaxies 
are "redder" on average (as can also be seen on Fig. 8, bot- 
tom right) . De Jong's data support this picture qualitatively. 
Also, the slope of the B-K vs. B-V relation {bottom right pan- 
nel in Fig. 12) depends little on A or Vc; it is slightly steeper 
than the observed one, but the overall comparison is fairly 
satisfactory. 

It is obviously difficult to assign an "age" to these galax- 
ies on the basis of such diagrams. Increasing slightly the ef- 
ficiency of star formation in Eq. (1) would lead to similar re- 



sults at ages lower by a few Gyr. Further constraints should 
be considered, like those analysed in the next sections. 

We wish to emphasize here the fact that our adopted 
senario (infall and star formation timescales longer for low 
mass galaxies) leads naturally to the trend of massive galax- 
ies being "redder" on average, as observed. This trend is a 
well-known problem in hierarchical models of galaxy forma- 
tion, where low mass galaxies are formed first and are found 
to be "redder" today than massive ones. In their recent work 
Somerville and Primack (1998) suggest that this problem 
may be overcome by including in the models dust extinction, 
which affects more the massive galaxies and makes them 
"redder" on average. However, in our models we find that 
dust extinction plays a minor role: it can modify slightly 
the situation (enhancing the already existing trend in the 
stellar population alone) but not reverse it completely. And 
is certainly not dust extinction that makes the Milky Way 
appear chemically and photometrically older than the Mag- 
ellanic Clouds. 
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Figure 13. Gas fraction vs. B-magnitude {left pannel), B-V colour {middle pannet) and central surface briglitness ^bo (Tight pannel). 
The grid of models represents 5 values for the spin parameter A {solid curves) and 5 values for the circular velocity {dashed curves), as 
indicated on the left pannel. Values of Vc run from 80 km/s {open triangle) to 360 km/s [filled triangle). Results of models are given at a 
galactic age of 13.5 Gyr. Data in all pannels are from McGaugh and DeBlok (1997). Filled symbols correspond to high surface brightness 
galaxies and open symbols to low surface brightness galaxies. Our grid of models clearly does not fit the lowest surface brightness galaxies 
(i.e. below ^so=23 mag arcsec"^), which probably require values of A >0.1. In all three pannels, the dashed diagonal line is a fit to the 
data given in McGaugh and DeBlok (1997): Og = 0.12 (Mi3+23) = -1.4 [(B-V)-0.95] = 0.12 (^tso-lS). 



4.4 Chemical Evolution: Gas fraction vs. 
Magnitude and Colour 

Diagnostics of galaxy evolution can be made in a variety of 
ways. Colour-colour diagrams are often used for that pur- 
pose, as discussed in the previous section. However, their 
usefulness is limited by the fact that they admit "degener- 
ate solutions", i.e. "red" stellar populations may result ei- 
ther from old age and low metallicity or young age and high 
metallicity. 

A key ingredient to galaxy evolution is the consumption 
of gas, turned into stars; as time goes on, the gas fraction (jg 
decreases in general (except in cases where infall exceeds the 
star formation rate). The gas fraction could then, in princi- 
ple, be used to probe the evolutionary status of galaxies. 

This probe has rarely been used in studies of galaxy 
evolution up to now, for several reasons. First, it is also af- 
fected by "degeneracy" : galaxies may have low gas fractions 
either because of quiescent star formation for a long time or 
because of recent enhanced star forming activity. Secondly, 
until recently, the observed range of values in ag was rather 
small (cTg <0.2) and provided little constraints to models. It 
seemed indeed that we live at a special epoch, when most 
galactic disks have almost exhausted their reservoir of gas 
(e.g. Kennicutt et al. 1994). 

In an important recent work McGaugh and deBlok 
(1997) compiled data from several surveys and derived 
various important quantities for galaxy evolution studies: 
colours, absolute magnitudes, scalelengths, central surface 
brightness, atomic and molecular gas amounts. On the basis 
of these data they were able to establish several important 
correlations between the gas fraction and other quantities, 
in particular B-V and B-I colours, absolute magnitudes and 
central surface brightness in the B- and I- band. On the other 
hand, they found no apparent correlation between gas frac- 
tion and disk scalelength. Another important finding of their 
work is that gas rich galaxies (i.e. with erg >0.4) do exist, 
but they are systematically low surface brightness galaxies. 

In Fig. 13 we compare our results to the data of Mc- 



Gaugh and deBlok (1997). The following points can be 
made: 

1) We find indeed a trend between gas fraction and 
absolute magnitude: smaller and less luminous galaxies are 
in general more gas rich, as can also be seen in Fig. 8 {bottom 
left pannel) . Our grid of models covers reasonnably well the 
range of observed values, except for the lowest luminosity 
galaxies (see also points 2 and 3 below). In particular, the 
slope of the o-g vs. M_b relation of our models is similar 
to the one in the fit to the data given by McGaugh and 
deBlok (1997). Moreover, among our models, closest to this 
fit are those with A=0.05 (and 0.04, if we interpolate to the 
next value of A=0.03), i.e. the value near the peak of the 
A-probability distribution that is obtained from numerical 
simulations. 

2) We find that the observed Og vs. B-V relation can 
also be explained in terms of galaxy's mass, i.e. more mas- 
sive galaxies are "redder" and have small gas fractions. In- 
deed, the slope of the (jg vs. B-V relation (-1.4) is again well 
reproduced by our models, and the A=0.05 models lead to 
results that match close the observations. However, although 
our grid of models reproduces well the observed dispersion 
in CTg, it does not cover the full range of B-V values. The 
lowest B-V values correspond to gas rich galaxies {og ~ 0.5- 
0.8) and could be explained by models with higher A values 
than those studies here (i.e. by very Low Surface Brightness 
galaxies). However, observed galaxies that are both "red" 
(B-V~l) and gas rich ((jg >0.3) are difficult to explain in 
the framework of our models. 

3) The situation is radically different for the ag vs. fj,Bo 
relation. As already seen in Fig. 8 {top right pa,nnel) we find 
no correlation in our models between fXBO and Vc, but we do 
find one between ^bo and A: more compact disks (lower A) 
have higher central surface brightness. It is this latter prop- 
erty that allows our models to reproduce, at least partially, 
the observed Ug vs. ^bo relation {right pannel of Fig. 13). 
It can be seen that for a given A, similar values of /ibo are 
obtained for all galaxy masses (or Vc), but when different 
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values of A are considered, a large part of the observations in 
the ag-fiBO plane is covered by our grid of models. Galaxies 
with high fiBO (20-21 mag arcsec"^) and low ag correspond 
to models with A ~0.02. As in point 2 above, however, we are 
unable to reproduce galaxies with both high hbo and ag > 
0.3. On the other hand, our models do not reproduce very 
low surface brightness disks, with fiBO >23 mag arcsec"^; 
as explained in the previous paragraph, this is a rather mi- 
nor problem, since such disks could result from models with 
A >0.1, not considered here. 

We should like to stress that it is the first time (to our 
knowledge) that a comparison of this kind of data is made 
to fully self-consistent models of chemical and spectropho- 
tometric evolution of galaxies. Although the agreement be- 
tween theory and observations is by no means perfect, it 
nevertheless suggests that the basic picture may be correct. 
The main result is that the ag vs. M_g and ag vs. B-V rela- 
tions can be explained in terms of galactic mass (or, equiva- 
lently, Vc), but the ag vs fj.BO relation can only be (partially) 
explained in terms of the spin parameter A. 

It is interesting to notice that in their work McGaugh 
and deBlok (1997) have speculated on these relations, on 
the basis of older, one-zone models of galaxy evolution 
(by Guiderdoni and Rocca-Volmerange 1987). In particu- 
lar, they anticipated that the galaxy's baryonic mass may 
be the driving factor in explaining the ag vs. Ms relation, 
as we find in our models. However, they suggested that the 
ag vs. fiBO relation could be explained if low surface bright- 
ness galaxies are on the average 3-5 Gyr younger than the 
high surface brightness ones. Instead, we interpret the ob- 
served correlation in terms of A, not of age differences. We 
are aware though that age is a third parameter that affects 
the results and may render the picture more complex than 
scetched above. We prefer, however, to keep a minimum of 
simplicity in our models and assume that all disks started 
forming stars at least 10 Gyr ago, unless there is compelling 
evidence for the contrary. 

4.5 Chemical evolution: Metallicity vs. 
integrated and local properties 

Metallicity is another probe of the evolutionary status of 
a galactic disk. Both its absolute value and its radial pro- 
file can be used as a diagnostic tool to distinguish between 
different models of galactic evolution. 

Zaritsky et al. (1994) provided a sample of 39 disk galax- 
ies for which oxygen abundances have been measured in at 
least 5 HII regions. They determined radial abundance pro- 
files and deduced the oxygen abundance at three caracter- 
istic radii: an absolute physical radius (at 3 kpc from the 
center) and two dimensionless radii, namely 0.4 po (po be- 
ing the isophotal radius) and 0.8 ps {ps being the disk expo- 
nential scale-length, i.e. Rd in our notation). Zaritsky et al. 
(1994) found a strong correlation between the abundances 
at each of these radii and the disk circular velocity Vc, and 
a slightly weaker correlation between abundances and abso- 
lute magnitude Mb. They also found a correlation between 
abundances and Hubble type, but the fact that mass (i.e. 
Vc) and Hubble type are correlated with each other among 
spirals makes difhcult to evaluate the role of the latter in 
shaping these correlations. 

In Fig. 14 we plot our results at the same caracteristic 




Figure 14. Oxygen abundances at various galactocentric dis- 
tances: at 3 kpc [upper pannels), 0.4 po [middle pannles), and 0.8 
PS [lower pannels). They are plotted as a function of Vc [left) or, 
equivalently Mb [right). Our results [curves) are parametrised by 
the corresponding A values. Data are from Zaritsky et al. (1994). 



radii as Zaritsky et al. (at 3 kpc, at 0.4 po and at 0.8 Rd, 
from top to bottom, respectively) as a function of circular 
velocity Vc {left pannels) and of Ms [right pannels). 

It can be seen that a correlation is always found in our 
models between abundances [Z) and Vc or Ms, stronger in 
the case of Z[3 kpc) than in the other two cases. Only in the 
A=0.01 models (unrealistic for disks, as stressed already), 
Z[0.8 Rd) is not correlated with either Mb or Vc- Except 
for that, the grid of our models reproduces fairly well the 
data, both concerning absolute values and slopes of the cor- 
relations. Again, it is the first time (to our knowledge) that 
a detailed comparison of this type of data is made to fully 
self-consistent models of galactic chemical and photometric 
evolution. We think that the results of this comparison are 
quite satisfactory. 

Zaritsky et al. (1994) point out that the trend of Z 
vs. Ms observed in their galactic disks is simply the con- 
tinuation of the abundance-luminosity relation noticed by 
Garnett and Shields (1987) over a wide range of magnitudes 
between spirals and irregulars. The relation spans over 10 
magnitudes in Ms and over a factor of 100 in metallicity. 
If galactic mass is always driving this relation (as in our 
models), then little room is left for galactic winds to play 
an important role; indeed, it is often claimed that low mass 
galaxies suffer loss of metals and/or gass through galactic 
winds, but in that case a deviation from the correlation ob- 
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Figure 15. Oxygen abundances as a function of radius (in units of the effective radius R^jjj, left pannel), local surface density (middle 
pannel), and local gas fraction (right pannel). Our model results (open sumbols) are given at an age of 13.5 Gyr and correspond to 
different galactocentric distances; they are parametrised by the corresponding A value, as indicated on the right pannel. Data (filled 
symbols) are from Shields et al. (1991) for HII regions in field spirals. 



served in high mass galaxies (presumably unaffected by mass 
loss) would be expected. 

It should be noticed that Zaritsky et al. (1994) find 
no correlation between metallicity and gas fraction ag in a 
limited sub-sample (14 galaxies) of their data. As shown in 
Fig. 8, we do find such a correlation in our results (higher 
Z for lower ag). However, this correlation is stronger for 
high ag, while at low ag (below ~0.2) metallicity tends to 
"saturate". Zaritsky et al. (1994) admit that their negative 
result may be due to the fact that their sub-sample is limited 
to relatively low gas fractions. It is interesting to notice that 
our results show also that metallicity tends to saturate at a 
value of log(0/H)-|-12 ~ 9.3 for the most massive disks, i.e. 
those with the lowest gas fractions according to Fig. 8. 

The results of Zaritsky et al. (1994) concern abundances 
as a function of global galactic properties. Abundances in spi- 
ral galaxies as a function of local properties were studied by 
Shields et al. (1991), who searched for environmental effects 
between field and cluster galaxies (those in Virgo). They 
compiled data on abundances in HII regions and atomic and 
molecular gas surface densities as a function of galactocen- 
tric radius in five field galaxies. We concentrate on those 
field galaxies only, that presumably evolved "passively" (i.e. 
with no environmental effects) and correspond better to our 
models. 

The data of Shields et al. (1991) are plotted in Fig. 15 
and compared to our results. Oxygen abundances are plot- 
ted as a function of radius R (normalised to effective radius 
Rsff), of total surface density E and of gas fraction ag. Our 
results are shown at an age of 13.5 Gyr, parametrised by the 
values of spin parameter A. It can be seen that: 

1) There is a good anti-correlation in general between 
metallicity and R/REff,'S and ag [left, middle and right 
pannels, respectively). The scatter in the observed relation 
is larger in the case of Z vs. R/Rsff and smaller in the case 
of Z vs. E. 

2) Our models reproduce well, in general, the observed 
trends. However, the scatter in the observed Z vs. R/Rsff 
relation is only reproduced when A=0.01 models are con- 



sidered: they are the only ones producing high metallici- 
ties at large R/Rejj- But such A values produce unrealistic 
disks, as stressed several times in this work ; and if we ne- 
glect them, we are unable to reproduce the observed scat- 
ter. The solution to that difficulty may reside in the fact 
that the large abundances at large R/Rsff in the Shields et 
al. (1991) data come from only one galaxy, N2903. If this, 
perhaps "peculiar" , galaxy is removed from the sample, the 
rest of the data present a tighter correlation that is readily 
explained by our models with A >0.02. 

3) The metallicity vs. surface density relation presents 
the tightest correlation, both observationally and in our 
models. We find an excellent agreement between theory and 
observations in that case. 

4) The metallicity vs. gas fraction relation is the most 
intriguing. Our model results are systematically ~0.2-0.3 
dex higher than the observations. One might think that the 
adopted yields are too high by ~50%; such an uncertainty in 
stellar yields is indeed conceivable (see Prantzos 1998). How- 
ever, in that case our results as a function of R/ Re / / and E 
should also be modified: lowering the theoretically predicted 
O/H by 0.2-0.3 dex would spoil the excellent agreement with 
the data. Another possibility is that gas fractions have been 
systematically underestimated in Shields et al. (1991): mov- 
ing the data points to the right (in the right pannel of Fig. 
15) could alleviate the discrepancy. Also, a combination of 
those two factors (i.e lowering the theoretical metallicity and 
increasing the observed gas fraction) would help. We wish, 
however, to point out an interesting parallel with the sit- 
uation in our own Galaxy. It has been noticed on several 
occasions (e.g. Pagel 1997) that solar CNO abundances are 
higher by 30-50% than those in nearby young stars and the 
local interstellar medium (ISM). Standard galactic chemi- 
cal evolution models are usually required to reproduce solar 
abundances 4.5 Gyr ago and they always find present day 
CNO abundances larger by ~30-50% than those observed in 
the local ISM (while reproducing correctly the present-day 
local gas fraction, see e.g. Prantzos and Aubert 1995). Our 
models being calibrated on the Milky Way, it is interesting 
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Figure 16. Spectra of our model disks at an age of T=13.5 Gyr, compared to observed spectra from the Catalogue of Kennicutt 
(1992). We have selected a few representative values of A (0.01, 0.03, 0.07 from right to left) and of Vc (80, 220, 370 km/s, from top 
to bottom) and plotted the corresponding integrated galactic spectra at T=13.5 Gyr {solid curves). For each of these cases, we selected 
in Kennicutt's catalogue spectra of galaxies with measured rotational velocities (Vmax in each pannel, following the galaxy's name) as 
close as possible to the model values; they are plotted in dotted curves. All spectra are normalised to the flux at 550 nm. Notice that our 
model spectra do not include emission lines. 



to see that they overproduce also present day oxygen abun- 
dances in field spirals, and by similar amounts! This consti- 
tutes one more indication that solar abundances may not be 
representative of the local ISM 4.5 Gyr ago, as pointed out 
in several works (Prantzos et al. 1996, Clayton 1997). Con- 
firmation of this would have major implications for galactic 
chemical evolution studies. 

4.6 Spectral evolution: integrated spectra 

We calculated the integrated spectra of our galactic disk 
models, according to the procedure presented in BP99. Our 
aim is to see whether our models can match reasonably well 
observed spectra of spiral galaxies. 

In Fig. 16 we show our results for three values of A 
(0.01, 0.03 and 0.07, from left to right) and three values of 
Vc (80, 220 and 360 km/s, from top to bottom) at an age 
of 13.5 Gyr {solid curves). As cen be seen, the form of the 
spectrum depends on both A and Vc- Low values of A and 
high values of Vc lead to old stellar populations and to a 
severe depression of the spectrum "bluewards" of 400 nm, 
due to an efficient early star formation (see Fig. 4 and 8). 

We searched in Kennicutt's (1992) catalogue of galactic 



spectra galaxies with measured rotational velocities close to 
our Vc values. For each one of them we were able to find 
observed spectra that match reasonably well our theoretical 
ones {dotted curves in Fig. 16). This is the case, for instance, 
of the large spirals NGC2775 and NGC3147 (K„a^=302 and 
333 km/s, respectively, in the middle and righ bottom pannel 
of Fig. 15). 

We stress that the successful comparison of observed 
spectra to our models proves nothing about the quality of 
the latter. Indeed, there are several ways to obtain spectra 
corresponding to young, intermediate and old stellar popula- 
tions, and use of composite spectra alone is a poor diagnostic 
tool (especially if extinction has to be taken into account). 

However, the comparison of models to observations 
made in Fig. 16 (the first ever made that involves mea- 
sured rotational velocities), when combined to other probes 
of galactic evolution (all those discussed in the previous 
sections) constitutes a powerful tool of galaxy diagnostics. 
Indeed, it would be interesting to see, for instance, what 
are the values of gas fraction, chemical abundances, colours, 
magnitudes, disk scalelengths etc. for the observed galaxies 
of Fig. 16, and check whether they fit to the general scheme 
elaborated in this work. Also, it would be interesting to find 
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out what makes galaxies with similar Vmax (like NGC3471 
and NGC 4775) display so different spectra; our models sug- 
gest that spin parameter may be the answer. In any case, a 
detailed study of this kind requires as uniform a set of data 
as possible. Work is curently in progress in that direction 
(Boissier and Prantzos, in preparation). 



5 SUMMARY 

In this work, we present a comprehensive study of the chem- 
ical and spectrophotometric evolution of spiral galaxies, us- 
ing multi-zone models. With respect to previous studies of 
that kind, our work presents the following important fea- 
tures: 

- The fully self-consistent treatment of chemical and 
spectrophotometric evolution, with up-to-date, metallicity 
dependent, stellar input data. 

- The use of a realistic stellar IMF, i.e. measured in 
the solar neighborhood (Kroupa et al. 1993), instead of the 
often used Salpeter IMF in the whole mass range; also, the 
radial dependence of the SFR (proportional to the rotation 
frequency V{R)/R, with realistic V{R) profiles), based on 
theories of star formation induced by spiral waves in disks. 

- The calibration of the model on the Milky Way disk 
(BP99), fixing the star formation efficiency which is not a 
free parameter for the other disk models. 

- The use of simple "scaling laws" for the properties 
(sizes and surface densities) of other galactic disks, that are 
adopted as boundary conditions in our study; these scaling 
relations are based on currently popular models of galaxy 
formation, in the framework of Cold Dark Matter senarios 
and allow to describe disks in terms of two parameters: ro- 
tational velocity Vc and spin parameter A. 

The main results of our models may be summarised as 
follows: 

i) The final gas fraction ag, M/L ratio, metallicity and 
colours depend on both A and Vc- Disks of low A (compact) 
and/or high Vc (massive) are formed earlier (in less than 6 
Gyr) and are "redder" today than disks with other {X,Vc) 
values. 

ii) Some disk properties are found to be quasi- 
independent of Vc (like central surface brightness fio, in 
agreement with observations), while others are independent 
of A (like e.g. current star formation rate 'I'o, absolute mag- 
nitude Ms or colour-colour relation). 

iii) High A values (>0.07) lead to Low Surface Bright- 
ness galaxies (as already shown in previous works, e.g. 
Jimenez et al. 1998), while A=0.01 leads to unrealistically 
small and bright disks. 

iv) Due to the inside-out star formation scheme adopted 
here, profiles of surface brightness and colours fiatten in 
time. 

Comparison of our results to a large body of observa- 
tional data leads to the following conclusions: 

1) Disk size and central surface brightness are well re- 
produce by our models, implying that the adopted scaling 
laws produce "realistic" boundary conditions; only A=0.01 
values produce unrealistic results, ressembling more to 
galactic bulges than to disks. 

2) The Tully-Fischer relationship is "built-in" in this 
type of models (i.e. through the scaling relations, which con- 



stitute the boundary conditions of the problem), but realistic 
senarios for star formation lead to small but important vari- 
ations: in particular, our models predict a scatter increasing 
with decreacing circular velocity (a trend suggested by the 
data of Giovanelli et al. 1997) and a steeper slope in the I- 
band for low surface brightness galaxies. Also, we find that 
extinction by dust makes the TF relation flatter than in 
the case of the stellar population alone. Most importantly, 
our models naturally predict an increase in the slope of the 
TF relation with decreasing wavelength, resulting from the 
fact that less massive galaxies are chemically and photo- 
metrically "younger" than more massive ones; this trend is 
indeed observed (TuUy et al. 1998). 

3) Colour-colour and colour-magnitude relations are 
well reproduced with the adopted (Milky Way-type) star 
formation efficiency, provided the observed disks are more 
than ~10 Gyr old. More massive disks are "redder" on av- 
erage than low mass ones. 

4) The relation of gas fraction to magnitude Ms, colour 
B-V and central surface brightness hbo is well reproduced. 
The trends in the former two correlations are driven by 
galaxy's mass (Vc), while in the latter by the spin param- 
eter A. Some of the observed low surface brightness galax- 
ies require models with higher A values than those studied 
here. Observed high surface brightness galaxies with large 
gas fractions {ag >0.3) cannot be explained by our models. 

5) Metallicity is found to be well correlated to mass 
and luminosity, in fair agreement to observations. Also, it is 
correlated to local environment (radius, total surface den- 
sity), again in agreement with observations. Our models 
show higher metallicity than observed at a given mass frac- 
tion, reminding of an analogous situation in the Milky Way 
disk. This "problem" may represent another clue against the 
usual assumption that the solar abundances are typical of 
the local ISM 4.5 Gyr ago. 

6) Observed spectra of spirals with measured rotational 
velocities are well reproduced by our models, when compa- 
rable values of Vc are considered. The spin parameter A can 
help to explain differences between galaxies with the same 
rotational velocity. 

The agreement of our models to such a large body of ob- 
servational data is impressive, taking into account the small 
number of free parameters. Indeed, the only really "free" pa- 
rameter is the dependence of the infall timescale on galactic 
mass (Fig. 3), since all other ingredients have been fixed 
already by the calibration of the model to the Milky Way 
disk. Although it does not constitute a proof of the validity 
of the model, this success suggests nevertheless a coherent 
overall picture for the evolution of galactic disks, along the 
following lines: 

- Baryonic gaseous disks form in non-baryonic haloes, 
with properties given by simple scaling relations of the Cold 
Dark Matter senario. 

- Timescales for star formation in the disks depend 
mainly on their mass: massive disks form earlier the bulk 
of their stars and reach higher metallicities and B-V values 
than their low mass counterparts. This picture seems to con- 
tradict the idea that low mass haloes (and baryonic disks) 
form first. However, the timescale for assembling the bary- 
onic gas need not be the same as the star formation time 
scale: low mass galaxies may assemble relatively rapidly, but 
require a very long time to turn their gas into stars. This 
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point is crucial, because we found that it allows to explain 
simultaneously three important observed features: the vari- 
ation of the slope with wavelength band in the TF rela- 
tion (Sec. 4.2, Fig. 11), the fact that more massive galaxies 
are on average "redder" than low mass ones (Sec. 4.3) and 
the metallicity-luminosity relation of Zaritsky et al. (1994, 
Fig. 14). It is important to notice that, the assumption of 
mass-dependent galactic winds could help explain the third 
feature, but not the first one. 

- Star formation in disks takes place inside-out, produc- 
ing steep early gradients in metallicity and colour profiles, 
that flatten at late times. 

Several other predictions of this senario can be checked 
against observations. In a forthcoming paper (Boissier and 
Prantzos, in preparation) we show that the resulting abun- 
dance gradients are in agreement with observations. More- 
over, the senario predicts little evolution of massive disks at 
late times (since most of the action takes place early on) 
and more important late evolution for low mass disks. The 
former point seems to be supported by observations of disk 
properties up to redshifts of ~1 (Lilly et al. 1998). 
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